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Bose-Einstein condensates of polar molecules: anisotropic interactions = anisotropic 

mass 
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So far the theory of Bose-Einstein condensates (BEC) of polar molecules was based on an ad hoc 
generalization of equations for spherical atoms. Here I adopt a rigorous pseudo-potential approach 
to low-energy dipolar interactions and derive a non-linear mean-field Schrodinger equation for a 
harmonically-trapped condensate. I show that, effectively, the dipolar interactions alter molecular 
mass. The resulting effective mass is anisotropic: to the leading order the mass is altered only for 
the motion along the polarizing field. For a typical BEC of spin-polarized magnetically-interacting 
alkali-metal atoms the effective atomic mass is reduced by 10% from it's bare value. For a BEC of 
polar molecules the mass may be reduced by a factor of a 1,000. 
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The ongoing quest for a Bose-Einstein condensate of 
polar molecules [l| is stimulated by a remarkable rich- 
ness of the quantum phenomena associated with their 
large dipole-dipole interactions. Inherently anisotropic, 
the dipolar interactions are crucial to quantum informa- 
tion processing Hand to realizing exotic states of matter 
(e.g., topological[3] and checker-board [4| phases) with ul- 
tracold molecules. 

Background — The intricacies of the many-body prob- 
lem are rooted in inter-particle interactions that lead 
to non-separable Hamiltonians. Introducing pseudo- 
potentials, i.e., effective interactions that are simpler 
than the original interactions, makes the problem more 
tractable. In physics of ultracold atoms, all seemingly 
disparate interactions can be well approximated by a con- 
tact pseudo-potential. Its strength is determined by s- 
wave scattering length a ss , which is computed by solving 
low-energy scattering problem with the original poten- 
tial. This contact interaction is central to arriving at 
the non-linear Schrodinger equation for the condensate 
wavefunction \P (r). The hallmark Gross-Pitaevskii equa- 
tion H (GPE) reads: 



2M 



C/(r)+^|v|/(r)| 2 ^( r ) =/i0 vl/( r ), (l) 



where [1q is the chemical potential, and U (r) is the 
external confinement potential. The non-linear term 
arises due to interparticlc interactions, g' = Anh jM a ss . 
Namely this non-linearity yields a wealth of non-trivial 
effects, such as solitons, in physics of the condensates 6]. 
Why do the dipolar interactions require going beyond 
the conventional approximation? Compared to the l/rf 2 
dependence of the conventional isotropic van der Waals 
interactions on interparticle separations, the dipolar in- 
teractions are both long-range, l/rf 2 , an d anisotropic. 



This seemingly innocuous power-law variation crucially 
modifies low-energy collision process that underlies the 
pseudo-potential formalism. Contributions of partial 
waves beyond the s-wave no longer "freeze out" and the 
scattering is characterized by an infinite number of scat- 
tering lengths. Due to the anisotropy, molecules exert 
torques on each other and various spherical waves (e.g., 
s and d) become coupled. These couplings result in ad- 
ditional "anisotropic" scattering lengths [7|. 

Following Yi and You@] (YY), the r apidly growing lit- 
erature on dipolar BECs, see e.g. HSISIlJlSliaEl 
LL5[ , is based on an effective interaction that is represented 
as a sum of the contact pseudo-potential and the classi- 
cal dipole-dipole interaction. This ad hoc approximation 
has a shortcoming of being valid only in a perturbative 
(Born) regime. Another, both aesthetic and practical 
shortcoming, is that in the YY approximation the GPE 
becomes a non-linear integro- differential eq\idXio\i\&\ that 
lacks the appealing minimalism of Eq. (fT]). By contrast, 
here, starting from the rigorous quantum-mechanical de- 
scription of the dipolar collision process Q , we attempt 
to overcome both shortcomings: (i) the employed pseudo- 
potential involves scattering parameters that may be 
tuned all the way through the resonances and (ii) the 
dipolar GPE derived here has a simple mathematical 
structure. 

Dipolar pseudopotential — In a typical dipolar BEC 
setup, a molecular gas forms a cloud in an external 
harmonic trapping potential. Orientation of molecular 
dipoles D is fixed by applying a polarizing E-field (other- 
wise, molecular rotations would average dipole moments 
to zero). Then as r\i — > oo, the molecular interactions 
acquire dipolar character, 
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Here d\i is the angle between collision axis Vn and the 
polarizing field. The collision process is also determined 
by the short-range part of the potential: as molecules ap- 
proach each other, the electronic clouds start to overlap, 



and the interactions substantially depart from the dipo- 
lar form @. The YY approximation treats the short- 
and long-range parts of the full interaction on separate 
footings; we will incorporate both consistently. 

We need to describe a quantum dipolar collision pro- 
cess at ultralow temperatures. I assume that the polar- 
ized atoms or molecules follow a unique potential surface. 
For example, recently attained BEC of highly-magnetic 
chromium is comprised of spin-polarized atoms [1 €31 ] . The 
atoms are trapped in the lowest-energy Zeeman sublevel; 
transitions to the upper-energy levels are forbidden ener- 
getically. To quantify the scattering, one has to solve a 
multi-channel problem. The relevant scattering parame- 
ters are the following limits of the K-matrix characteriz- 
ing couplings between I and £' partial waves pj], 



lim K,( m - 



,/k. 



(3) 



where fik is the relative momentum of the colliding pair. 
K-matrix essentially governs the asymptotic form of the 
scattering wavefunction for large interparticle separa- 
tions. Long-range, l/rj 2 , character of the dipolar in- 
teraction ensures that the above limits are finite. The 
quantities ag m -e> m ' have a dimension of length and will 
be referred to as scattering lengths. 

For illustration, consider scattering lengths for a pure 
dipolar interaction (i.e., assuming the validity of Eq. ((2J) 
for all 7-12) in the Born approximation. The formalism 
is described in Ref. Q and it is based on a system of 
coupled radial equations for individual partial waves. We 
find that in the Born approximation both the diagonal 
and off-diagonal scattering lengths fall off as t~ 2 with 
increasing £. In the following I assume that the dominant 
effects are due to s — s and s — d scattering lengths, a ss 
and a s d- The former is mainly determined by the short- 
range part of the potential and to the latter by the dipolar 
coupling. In the Born approximation [7[ 



Born 
a sd 



1/(6^5) MD 2 /Ti 2 



(4) 



Values of af£ ln for molecules of present experimental in- 
terest are listed in Table [TJ The off-diagonal scatter- 
ing length is strongly suppressed (by ~ (1/137) 2 ) for 
magnetically-interacting atoms: it is -0.01 nm for Cs and 
-0.2 nm for Cr 0. 



OH (X tI 3/2 ) 
RbCs (X X S) 
KRb (X X E) 
NH (A 3 S~) 



D, Debye af d orn , nm 



1.65 
1.2 
0.59 
1.38 



-52 
-350 
-48 
-32 



TABLE I: Anisotropic scattering length a a d in the Born ap- 
proximation for molecules of current experimental interest [1(. 



Pseudopotential — The low-energy pseudopotential for 
anisotropic scattering was introduced in Ref. [7| ■ An ear- 
lier variational BEC study with this pseudopotential (al- 
beit its erroneous version) was carried out in Ref. [17| . 



Recently, there was a study of the validity of the pseu- 
dopotential approach [l8| • These authors found that the 
pseudopotential description remains accurate as long as 
a sc i is smaller than the characteristic length of the trap- 
ping potential. 

For the goals of this paper we simplify the rigorous 
pseudopotential [7] . More details of the discussion pre- 
sented below can be found in Ref. |19(. We assume 
that the global BEC properties can be described by 
well-behaved wavefunctions. In this case we may op- 
erate in terms of the momentum-space representation. 
The matrix element of the pseudopotential between two 
plane waves (r|k) = (2tt) ' c' kr is given by v (k. k') of 
Ref. Q. For our case of the dipolar interactions trun- 
cated at s — s and s — d couplings it reads €> (k, k') = 



1 /il 

2tt' 2 M 



(a ss - a sd T (k, k')) with 



^(k,k') - Vb{p 2 (cos0 fe ) + {k/k') 2 P 2 (cos**)} , 

where 9k and 9k' are angles between the polarizing field 
and k and k'. 

Now, under simplifying assumption of harmonic trap- 
ping, I transform the momentum-space expression back 
into the coordinate space. We write for a matrix element 
of the pseudopotential (cf. Ref. [2(| for Rydberg atoms) 

(ip\V ps \ip) = (27i-r 3 f dkdk'drdr' x 

V>* (r') e lk '- r 'v (k', k) V (r) e^ r . (5) 



Only certain values of |k| and |k'| contribute to this 
integral. Experimentally, the collisions occur in the 
presence of harmonic trapping potential, U (r) = 
\M {uj 2 x 2 + u> 2 y 2 + u 2 z 2 ). For two harmonically- 
confined particles the center-of-mass and relative motions 
decouple and the Hamiltonian for the relative motion 
reads H r = p\ 2 / (2/x) + (&) U (|r 12 |) + V(r 12 ), where 
V (ri 2 ) is the full interaction potential between the par- 
ticles and /1 = M/2. In the stationary problem we solve 
the eigenvalue equation H r ip (ri 2 ) = E r if) (ri 2 ), E r being 
the energy of the relative motion. Ref. [21| presented a 
numerical comparison of solutions of this equation with 
the full molecular potential and its pseudopotential rep- 
resentation (for isotropic s-wave scattering). They ob- 
served that for velocity-dependent potentials the relevant 
collision momentum k c is h 2 k 2 /(2fi) — E r . Indeed, the 
collision process occurs at |ri2 1 much smaller than the 
harmonic length. In this region [/(|r 12 |) ~ and the 
kinetic energy is E r . This leads to k 2 = 3/2 Moj/h, with 
Cu = '^ /i LUi/3 being the average of the three trapping fre- 
quencies. 

Returning to the evaluation of the integral (JSJl , we see 
that the relevant contributions are accumulated at |k| = 
|k'| = k c . Representing cos 6^ = k z /k, and using the 



properties of the Dirac (5-function, we arrive at 



V ps (na) 



47r — S (ri 2 ) x 



(6) 
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a s d 
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This is a contact interaction that depends on the tra- 
ditional s-wave scattering length a ss and "anisotropic" 
scattering length a sd . The latter is due to the dipolar 
coupling of s and d partial waves. Both scattering lengths 
are to be determined from multi-channel scattering calcu- 
lations with the full interaction potential. One recognizes 
the conventional contact interaction AttH /M 5 (vi2)a as 
for the s-wave scattering; the remaining part, with a sd , 
is due to the dipolar interactions. In this term the 

9? = d 2 /dz 2 2 acts on the bra and 9% operates on 

the ket. The dipolar contribution to Vp S breaks into 
two parts: isotropic and anisotropic (derivative) terms. 
The isotropic dipolar contribution merely renormalizcs 
the traditional spherically-symmetric pseudopotential. It 
is the derivative term that defines the novel physics of the 
dipolar BEC. The anisotropy is imposed by the polariz- 
ing field (the derivatives are taken along the field). The 
long-range character of the dipolar interactions is mani- 
fested through k 2 which characterizes the entire trapping 
potential. 

The low-energy dipolar interactions can be controlled 
by external fields. As in the case of Feshbach-resonance 
tunability of a ss , the anisotropic length a sd may ex- 
hibit resonances. There is an evidence for such a 
field-dependent resonance in dipolar collisions of E-ficld- 
polarized Rb atoms [22j. This opens an intriguing pos- 
sibility of resonantly controlling and enhancing dipolar 
interactions in BECs. The described resonant scatter- 
ing mechanism is to be distinguished from a recent pro- 
posal [23| on controlling the strength of the dipolar in- 
teraction by modulating the dipoles by oscillating fields. 

Mean-field approximation — With the derived pseudo- 
potential, I proceed to analyzing properties of the dipo- 
lar BEC. In the mean- field approximation, all identi- 
cal bosons occupy the same single-particle wavefunction 
^f(r). Under this assumption, we arrive at the total en- 
ergy of the condensate (\E r is normalized to the total num- 
ber of particles N) 



Em 



dv 
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|W| 2 + [/(r)|*| 2 - 



\gort-J±\*\*] 



(7) 



where the isotropic and anisotropic coupling parameters 
are defined as 

.go = 4vr?i 2 /Af (a ss + V5a sd ), (8) 

g d = 2Trh 3 /(M 2 uj)V5a sd . (9) 



The dipolar effects are governed by the anisotropic length 
a sd . Dipolar interactions modify the traditional GP term 



ho 1*1 



.,.;,, ... i and appear in the newly-introduced derivative 
term, -g d (d z \$\ 2 ) 2 . 




aspect ratio k 



FIG. 1: Comparison of the dipolar interaction energy com- 
puted in the traditional YY approximation (solid line) and 
within the present pseudo-potential approach (dashed line). 
The comparison is given as a function of the aspect ratio of 
the cloud. 

In the present analysis, we truncated the expansion 
over partial waves at the s — s and s — d scattering 
lengths. To estimate the truncation error, we turn to the 
regime where the Born approximation is valid. Then we 
use Eq.Q and compare the results with those computed 
within the conventional YY approximation. In partic- 
ular, consider an important expectation value (Vdd) of 
the dipolar interaction for an axially-symmetric ground 
state of the harmonic trapping potential (u> x — w s = 
uj±, aspect ratio k = (Wz/wj,) 1 ' 2 , polarizing field is 
along the axis of symmetry). We may parameterize 
{Vdd) = N?kVx(afS Tn /ax) X v(k). In the YY approxi- 
mation v YY (k) — —6k/(\/2tt) b (k), where function b (k) 
is given in Ref. [8]. In the present pseudo-potential ap- 
proach v ps (ft) = \/^( K ~ r^+2l ) • ^ comparison of 

the two reduced dipolar energies v ps and v YY as func- 
tions of the aspect ratio is presented in FigfTJ We find a 
good agreement between the two curves. Both energies 
vanish for spherically-symmetric [k = 1) traps. 

Effective anisotropic mass — For real- valued 5", we 
combine the kinetic-energy and the derivative terms in 
Eq. ||7J). The result suggests introducing effective mass 
along the polarizing field, 



Mf z (r) = M/(l - 8g d M n (r) /tf) , 



(10) 



where the number density n(r) = \^{r)\ 2 . The mass 
remains "bare" (M) for the motion perpendicular to the 
polarizing field. In the Born approximation, the relative 
change in the mass is 



sm: 



cir 



M 



8nn{r)D 2 
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i.e., the mass is reduced from it's bare value by the ra- 
tio of the characteristic dipole interaction n(r)D 2 to the 
trapping energy. As an illustration, consider a BEC of 
alkali- metal atoms (magnetic moment 1 (J-b) with n rs 
10 14 1/cm 3 in a trap of u) = 2n x 100 Hz. Wc find that 
the mass is reduced by 10%. For a BEC of molecules the 
effect is even more pronounced. Here the effective mass 
is about a 1,000 times smaller than the bare mass. 
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mass 




FIG. 2: Effective anisotropic mass M° 2 (r) as a function 
of position for an "anomalous" (o a d > 0) BEC in a har- 
monic trapping potential. An unperperturbed Gaussian den- 
sity profile of cylindrical symmetry about the polarizing filed 
is assumed, n(r) = n max exp{ — (p/L±) 2 — (z/L z ) 2 }. The 
peak density is choosen to exceed the critical value, n max > 
h 2 /(8Mgd), and the effective mass M zz is negative in the 
center of the cloud. As we move to the outer regions of the 
cloud, n decreases, M zz goes through a singularity and be- 
comes positive. For even smaller densities, M zz approaches 
the bare mass M. 

The effective mass depends on the local density. Es- 
pecially interesting is the case when a sc i is tuned to 
positive values: here the effective mass may become 
negative (see Fig. [2]). Here the condensate is unsta- 
ble in the region of negative mass: the energy may be 
lowered by accumulating density oscillations of increas- 



ingly shorter wavelengths. By requiring that at the 
peak density M°f > 0, we find that the BEC is sta- 
ble for a number of mo l ecules belo w the critical number 
N < N cr i t ~ l/8y/Tr/5y/h/(AI u> z )/a,sd- Such an instabil- 
ity is not present in the traditionally-considered dipolar 
gases, as af£ ln < 0. We observe that the instability in 
this case is likely to be related to the well-known instabil- 
ity due to attractive isotropic interactions [6] (when the 
effective scattering length a ss + V5a s d < 0) . 

Finally, we minimize Eq. ([7]) with respect to \&, and 
arrive at the non-linear Schrodinger equation 



2M 



A + J7(r)+flo|#(r)f *(r) 



8P 

dz 



9d {d^^ {r)l J * (0 = W* (r) . (11) 
This equation subsumes the traditional GPE |T]) when 
o-sd = 0, i.e., for bosons interacting via spherically- 
symmetric forces. Dipolar length a sc [ modifies the 
isotropic term and it governs the derivative term. By con- 
trast to the so far employed YY approximation, Eq. fTTj) 
is more concise and remains valid even in a vicinity of 
scattering resonances. 

To conclude, here I developed a new framework for 
analyzing dipolar BECs. The principal results are: the 
dipolar pseudopotential ©, the energy functional 0, 
and the mean-field equation (jll[) . The results may be 
interpreted in terms of the effective anisotropic mass: I 
showed that the interactions between dipoles alter mass 
for a motion along the polarizing field. For a typi- 
cal BEC of spin-polarized magnetically- interacting alkali- 
metal atoms the effective mass is reduced by 10% from 
it's bare value. For a BEC of polarized heteronuclear 
molecules the mass may be reduced by a factor of a 1,000. 
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